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ABSTRACT 

We study the survival of sub-structures (clumps) within larger self-gravitating dark matter 
haloes. Building on scaling relations obtained from A^-body calculations of violent relaxation, 
we argue that the tidal field of galaxies and haloes can only destroy sub-structures if spherical 
symmetry is imposed at formation. We explore other mechanisms that may tailor the number of 
halo sub-structures during the course of virialisation. Unless the larger halo is built up from a 
few large clumps, we find that clump-clump encounters are unlikely to homogenise the halo on 
a dynamical timescale. Phase-mixing would proceed faster in the inner parts and allow for the 
secular evolution of a stellar disc. 
Introduction 



High-resolution simulations of structure forma- 
tion in an ACDM cosmogony have revealed a 
large number of substructures within dark mat- 
ter haloes (e.g. Ghigna et al. 2000; Moore 2001). 
The mass function of these dark clumps (Moore et 
al. 1999a; de Lucia et al. 2004) 
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in the mass range lO'^ to 10^^ solar masses. These 
orbiting self-bound dark matter satellites would 
drag along baryonic matter to a non-negligible 
fraction of their mass. Consequently galaxies 
should harbour a large number of dwarf galaxies. 
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when only a handful are found (Kauffmann et al. 
1993; Moore et al. 1999a; see Binney & Merryfield 
1999). 

The situation is made worse from a dynamical 
standpoint. Massive dark clumps would perturb 
the vertical structure of galactic discs through 
tidal heating, when their thin structure suggests 
the immediate neighbourhood of discs is devoid of 
such perturbations (Toth & Ostriker 1992; Moore 
et al. 1999a). The large-mass end of the 
clump mass distribution function is robust 
in view of Simulation checks find the large- 
mass end of the clump mass distribution 
function robust to numerical resolution is- 
sues (Moore et al. 1999b; Gao et al. 2004; 
Power et al. 2003). Therefore, the orbital dis- 
tribution of clumps in phase space must allow for 
long periods of unperturbed evolution by the 
disc, as emphasised by Navarro (2002) and Font 
et al. (2003). A more severe problem is the nar- 
rowness of streams of stars associated with dwarf 
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galaxies tidal debris (Ibata ct al. 2001; Johnston 
at al. 2001), which suggests a smooth background 
halo potential in order to preserve the cohesion of 
the stream, in direct conflict with computer sim- 
ulations of structure formation. Observational 
evidence drawn from solar neighbourhood 
kinematics points to a fine-meshing of low- 
mass streams to account for coherent mo- 
tion in and out of the galactic plane (Helmi 
et al. 2002; Gould 2003). To settle this issue 
requires both a precise map of the mor- 
phology of haloes and convergence in the mass 
distribution function of clumps with a resolution 
down to the sub-dwarf mass range, still a 
challenge to present-day computer models (see 
e.g. Power ct al. 2003). The morphology of 
haloes has been discussed at length by sev- 
eral authors (e.g., Moore 2001; Ghigna et 
al. 2000; Fukushige & Makino 2003; Power 
et al. 2003; Navarro et al. 2004). Several 
contributions in Natarajan (2002) give a 
broad overview of the formation processes 
and equilibrium properties of haloes de- 
rived from numerical simulations and pos- 
sible observational tests for their detection. 

Particle-based calculations of galaxy formation 
proceed with a number routinely approaching a 
few xlO'' particles for a whole simulation. The 
mass of individual particles that make up dark 
haloes in these simulations ~ IO^Mq is still large 
compared with a mass spectrum of dark matter 
that may yet include a population of brown dwarfs 
(i.e., stellar masses). Furthermore, of order 
~ 10^ mass elements must participate in the 
formation of a self-gravitating body to re- 
solve the growth of potential energy ade- 
quately (Boily et al. 2002; Roy & Perez 
2004), or a mass resolution of 10^/10^ ~ 
10'^ Mg to account for dwarf-size structures. 
This leaves a gap in mass resolution of some three 
orders of magnitude with that achieved by present- 
day simulations. In this context, we need estab- 
lish when the statistics of dark matter clumps de- 
rived from N-body computations may be scaled 
up to actual galactic systems. A crude picture of 
halo formation divides the process in two stages: 
one of rapid collapse on a free-fall time scale, fol- 
lowed by a second, longer period of sporadic accre- 
tion (Bower 1991; Lacey & Cole 1994; Zhao et al. 



2003a, b). The first stage involves structures 
spread over a narrow range of masses, while 
the second stage sees low-rate accretion by 
a dominating, central body. The growth of 
potential energy during the free-fall stage of for- 
mation sets the maximum phase-space density of 
the distribution function by efficiently, if incom- 
pletely, redistributing binding energy between par- 
ticles (Lynden-Bell 1967; van Albada 1982). 

The growth of potential energy during collapse 
of self-gravitating systems is a function of the ini- 
tial morphology and the number of particles (or 
quivalently, mass resolution) used in the compu- 
tation (e.g. Boily et al. 2002 and references 
therein). A direct consequence of this is that 
the maximum tidal field oc V^^(r) experienced 
by galactic satellites (dark or baryonic) as 
they fly through the system varies with the 
morphology of the underlying halo, i.e., its for- 
mation history. For instance, at fixed resolution, 
the strength of tidal fields reaches a higher max- 
imum for an initially spherical distribution than 
for an axi-symmetric or triaxial one. The simpli- 
fication of symmetric distributions allows to scale 
up the results with particle number exactly (see 
§2 below). Would the tidal heating experienced 
by dark clumps during infall, scaled up to actual 
galactic halo particle numbers, be sufficient to un- 
bind them? In this short contribution, we apply 
the results violent relaxation studies to the 
tidal heating of galactic substructures in a hier- 
archical Einstein-deSitter universe. We show that 
the tidal field developing during infall may yet be 
sufficient to erase sub-structures smaller than a 
critical linear size Ic if the particle number N is 
sufficiently large and the galactic halo shows axial- 
or spherical symmetry during infall. 

2. Scaling of tidal fields 

Aarseth, Lin and Papaloizou (1988) (hereafter 
ALP-l-88) have shown that the growth rate of 
global modes of fragmentation during violent re- 
laxation is such that the minimal radius achieved 
by a spherical distribution scales with the simula- 
tion particle number N as 
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where -R(O) is the initial system size enclosing N 
identical mass elements and R{tff) the radial size 
at the free-fall time tg defined by 



tff 



(3) 



32G <p(a,0)> 

where < p{a, 0) >= AnM/Sa^ is the mean den- 
sity inside the particle's initial radius, a. A re- 
cent study extended the result ^ empirically 
with N-body calculations to axisymmetric (cylin- 
drical) and triaxial initial configurations (Boily et 
al. 2002). ft was found in these cases that C scales 
with particle number as 

^ I N^/^ (axisymmetry) 
I constant (triaxial). 

The factor C achieved by systems starting from 
triaxial distributions remains < 40 for ~ 86% of 
the parameter space of axial ratios a : b : c, how- 
ever the maximum achieved by individual realisa- 
tions is highly sensitive to the initial axes ratios 
and may yet diverge. N-body simulations using 
direct-integration algorithms (ALP-l-88; Theis & 
Spurzem 1999; Boily, Clarke & Murray 1999) or 
FFT integrators (e.g., Boily et al. 2002) repro- 
duce these relations over several decades in par- 
ticle number, giving confidence that the results 
are well recovered independently of the numeri- 
cal method used. A consequence of ([2]) and ([4]) 
is that a spherical or axisymmetric distribution 
will contract ever more as the number N of parti- 
cles is increased. The morphological evolution of 
haloes and their substructures will not be well re- 
solved if at formation the particle number which 
take part in the relaxation phase is too low. Zhao 
et al. (2003a, b) show that the formation history 
runs through a rapid accretion phase, followed by 
slow growth of the outer halo. The scaling law 
([2]) and ([4]) obtained for violent relaxation would 
therefore apply best to the first, early phase of 
formation, but not the late one. 

3. Cosmological application 

The mean tidal field of a collapsing halo is eval- 
uated from the time-dependent second derivative 
of the gravitational potential 



F, = L 



X V^$|^ = -I, 



GM 



with r < the system radius R, and double- 
differencing with respect to r at fixed time yields 
a measure of the tidal field at R acting on a clump 
of size Ic- It is clear that the tidal force is unbound 
if during collapse R{tff) 0, as would occur e.g. 
for a large-N spherical distribution, cf. ([2]). 

Bound clumps would survive violent relaxation 
if their binding energy is higher than the tidal 
heating they suffer during that phase (van Albada 
1982; Tsuchiya 1998). As shown from (O, rem- 
nant structures will be severely disrupted if the 
maximal potential depth achieved is large. To 
progress further we need invoke a result for struc- 
ture formation in an expanding Universe. Present- 
day data support an asymptotically flat metric for 
the Universe, and hence the Einstein-deSitter cos- 
mogony remains attractive. 

3.1. Background tidal heating 

In an EdS universe the relation between a 
bound structure's mass and virial radius is (Kaiser 
1986; Padmanabhan 1993; Somerville & Primack 
1999) 



R cx A'P 
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with 7 known in terms of the power-spectrum of 
density fluctuations P{k) at wavenumber fc. The 
classic CDM power spectrum at the time r of 
galaxy formation accounts for the time-evolution 
of structure from a bottom-up point of view. 
The relation between P(fc; r) and the Zeldovich 
spectrum arising from post-inflation decoupling, 
P{k;ti), is (see Bardeen et al. 1986, Appendix 
G): 

a{ti 



P{k-T) 



a{t) 



T'{k)P{k;U) (7) 



where a is the cosmical expansion factor and T 
the tranfer function which is well approximated 
analytically by 



T{q) 



ln(l + 2.34g) 
2.34^ 



(1 + 3.9q + 259(7^ + 163q^ + 2027^" 

(8) 
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with q = \k\VO/{nh'^Mpc^'^), 9 = prei/lMp^ is 
the ratio of relativistic paticles to photon energy 
densities. On the largest scales the Zeldovich spec- 
trum P{k; t) cx P{k\ ti) cx fc is recovered from ([8|), 
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while on small scales P{k;T) cx k~^. If we fit 
the power spectrum locally to a power-law of in- 
dex Uk, P{k) oc the index shifts progressively 
from rifc = +1 to Uk = —3 as we explore smaller 
scales. The indices and 7 are linked through 
the mass-radius relation ([6]) 

7 = ("fc + 5)/6 . 

Note that 7 is related to the power-index i/ relating 
mean density and mass, p cx , by 1/ = 7 — 3 = 
{uk — 13)/6; the free-fall time oc l/y^ therefore 
scales with mass as 

iff (X M(13-»'=)/12. 

On the smallest scales rik — > —3 and te cx M"^/^, 
a steeper relation than on the largest scales when 
rife and tg cx A/. Therefore small clumps 

have fully virialised when the violent relaxation 
phase of the larger halo begins. 

Since there is no fixed scales of mass or radius 
in gravitational dynamics, all virialised structures 
obey the same relation ©. If we lump together 
all those of virial radius < r, of mass m chosen 
such that M/m = N ^ 1, we have 



R ( M 
- cx — 

r \ m 



(9) 



In a hierarchical universe, small structures form 
first and hence the N small clumps have virialised 
well before the large underlying halo, of total mass 
M . At constant mass the virial theorem provides a 
relation between equilibrium radius an initial size 
Ri for a self-gravitating system: 



Ri — IR 



(10) 



which applies equally to all structures. We need 
relate the tidal field © with ^ and (g]) to de- 
termine whether a structure of size survives 
the formation of the larger halo, of radius Ri. 
The energy transferred to a small clump by tidal 
forces during infall is adequately quantified by 
Spitzer's (1958) impulse approximation formula 
even for relatively slow encounters (Aguilar & 
White 1985). This gives confidence that it will 
hold in the present context, where velocities are 
high. Tsuchiya (1998) finds this to be correct in 
his study of relaxing Plummer distributions. The 
tidal energy gained AE by a substructure of mass 



m, radius r, and internal binding energy E may be 
evaluated for a single passage at velocity V across 
the background halo (or galaxy) potential at the 
time of collapse, when the tide is maximum. In 
the impulse approximation, this is given by 



AG'^lVPra 
iR{tsYV^ 



C^ — E. 

m 



(11) 



Note that Ijlip applies to the tidal field of the 
background potential, and does not account for 
individual clump encounters. These are discussed 
later. Substituting ([9]) in pT|) we get 



E 



(12) 



whence we deduce that AE ^ E ii the coefficient 
j\r-(n.fc-(-3)/2^3 rejj^aijis small. For spherical sys- 
tems the collapse factor C obeys © and therefore 



AE 



7V-("fc + l)/2 



(spherical distributions) 



while for axisymmetric or triaxial distributions we 
find from ^ 



AE j7V-("^-+2)/2 (spheroidal) , , 
(X < . (13) 

(triaxial). 



E 



JSf-ink+3)/2 



The consequences of these results in relation to the 
power-spectrum of density fluctuations is clear: 
for clumps orbiting in a collapsing spherical halo 
or galaxy, tidal heating will be ineffective provided 
rife > —1. If the underlying distribution is axi- 
symmetric, tidal heating will be ineffective when 
Uk > —2. However for triaxial initial conditions 
this will hold true if > —3. Since the power- 
spectrum of observed matter distribution is never 
steeper than — —3, we deduce that the bulk of 
substructures or clumps evolving in larger struc- 
tures, such as dark matter haloes, will survive the 
formation of triaxial larger structures, if gravity 
alone fixes their binding energy. 

3.2. Tidal heating due to other substruc- 
tures 

The above conclusion only concerns the re- 
sponse of clumps to the background tidal field. 
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We may consider the interaction between clumps 
themselves as they cross the dense system. To this 
end we consider the tidal heating by two equal- 
mass substructures during an encounter. Substi- 
tuting M ^ TO in pT|) and the radius R{ts) by the 
mean distance between clumps x ~ i?(tff)/iV^/'^, 
an effective impact parameter, we find 



AE = 



4G- 



2-^4/3 ry E. 

3M \RJ 



3(i?/iVl/3)4y2 



In the above we used V"^ « 2GM/R, with 
M = Nm the total system mass as before. The 
expression reduces to 

AE^-N^/^(-)^ E. (14) 
3 V i? / 

Clearly to achieve AE « E requires r ~ i? or 
iV 3> 1. we may simplify (jl4p by substituting for 
f using jo]). We then find 

3 

Thus for any appreciable number N the tidal heat- 
ing due to encounters between clumps will be sig- 
nificant if rife < —13/3 = ~4i. No regime of the 
power spectrum covers that range and hence en- 
counters between clumps never produce significant 
tidal heating. A direct consequence of this is 
that while haloes form at different redshifts 
and sample different regimes of the struc- 
ture power-spectrum, the mass distribution 
function of substructures should be robust 
against cut-offs or significant changes to 
its shape. It is not clear yet whether the 
scale-free nature of the clump mass func- 
tion ([T|) measured in N-body simulations 
(e.g. Ghigna et al 2000) can be extended 
to very small masses (see Gao et al. 2004). 

There is of course one situation when substruc- 
tures can heat-up one another through tidal forces, 
which is when f ~ -R(O) or i?(<ff) and iV ^ 1. In- 
deed when f matches the mean intcr-clump dis- 
tance « R{ts)/N^/^, we compute AE/E ^ 1, al- 
ways. However this situation is more appropri- 
ate to galactic mergers than the formation of halo 
through accretion of several sub-units, as would 
occur in any bottom-up calculation of galaxy for- 
mation. 



The survival of substructures is in part due to 
the large relative velocity V established under the 
mutual potential of all clumps. Survival during in- 
fall is no guarantee that the substructures would 
remain once the halo has achieved virial equilib- 
rium. To estimate more precisely the net rate of 
heating on one clump due to the background tidal 
field as it crosses the system is made difficult be- 
cause of the large changes in potential taking place 
during violent relaxation. This is best done with 
numerical N-body calculations tailored for this 
problem. The very high-resolution simulations 
performed by e.g. Ghigna et al. (2000) demon- 
strate the likely survival of most substructures in 
and around dark matter haloes post-virialisation, 
in support of the basic argument outlined here. 

4. Discussion and conclusion 

Galactic satellites survive the formation phase 
of triaxial haloes and galaxies and will be de- 
stroyed on long timescales as they orbit the host 
galaxy. This result is obtained both from large N- 
body simulations (see also e.g. Bullock et al. 2000; 
Ghigna et al. 2000) and from semi-analytic ar- 
guments (Moore et al. 1996; Taffoni et 
al. 2003) as well as the fiuid calculation pre- 
sented here, and is therefore robust. The initial 
conditions and subsequent evolution of N-body 
computer simulations still plague their interpreta- 
tion and application to observed galaxies in terms 
of simple estimates of satellite disruption times. 
For instance, Font et al. (2003) and Ardi et al. 
(2002) have questioned the rate of disc heating 
by in-falling dark satellites. These authors find 
that thin discs may yet remain stable despite a 
high count of bound dark matter clumps, pro- 
vided the clumps do not follow near-radial or- 
bits. The problem of disc heating would seem- 
ingly not occur if the inner region of the halo 
were isotropic in phase-space. We already noted 
that the destruction of dark satellites would be 
more effective in the deep potential of the inner 
halo. Recently Gao et al. (2004) have re- 
derived statistics of halo substructures in 
computer simulations and found them to be 
less concentrated than the host halo. This 
and the coherence of cold-stream satellite 
debris cold-stream debris could be inter- 
preted in the light of the present analysis 
as pointing to a near-spherical (and hence 
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destructive), early phase of halo formation. 
Several studies have argued for a more 
spherical morphology in the inner region 
of haloes (e.g., Blumenthal et al. 1986; Du- 
binsky 1994). The cooling of baryons at the 
heart of DM haloes would provide a mech- 
anism for this by locking the inner halo 
morphology to a rounder shape than ob- 
tained from strictly gravitational evolution 
(Dubinsky 1994; Prenk et al. 1996). The 
inner morphology of galactic haloes would 
not automatically be spherical if baryons 
have had time to cool and form discs before 
the halo assembly is comlete : when that 
is the case, the halo's inner morphology in 
equilibrium is even more sensitive to the 
formation history and depends for instance 
on the orientation of the discs as they merge 
(see Kazantzidis et al. 2004) 

Possible mechanisms that may disrupt dark 
clumps on the dwarf galaxy mass-scale and below 
include supernova blow outs through gas irradi- 
ation and expulsion (Efstathiou 1992; Somerville 
2002; Gnedin & Zhao 2002). The net effect of 
gas loss on gravitationally bound structures is un- 
likely to be effective if the gas mass fraction is 
small (Hills 1980; Boily & Kroupa 2002). Gnedin 
& Zhao (2002) have argued that the peaked den- 
sity profiles obtained from CDM numerical calcu- 
lations would resist rapid removal of the gas un- 
der any realistic circumstances. Thus unless the 
dwarf-size clumps contain a very large fraction 
of baryonic matter they will survive any degree 
of gas heating. Cote et al. (2002) presented a 
Monte Carlo simulation of chemical enrichment 
of galaxies through gas-evaporation from clumps 
of dark matter initially seeded with baryons (uni- 
form M/L ratio). They find that the mass func- 
tion of seeded clumps required to match their 
sample of galaxies (in terms of chemical gradi- 
ents and observable dwarf galaxy and star cluster 
populations) is similar to the mass function ([1]) 
obtained from large-N cosmological calculations. 
This would suggest that the high-mass end of the 
clump mass function survives the formation of the 
host dark halo to produce the observed popula- 
tion of dwarfs. It does not however suppress the 
number of dark clumps that may still be orbit- 
ing the halo. Another route to solving the over- 



abundance of dwarf galaxies is by preventing 
(bright) baryons from forming stars. This 
can be achieved either through background UV 
radiation (Efstathiou 1992; Somerville 2002), or 
by preventing a Toomre instability from devel- 
oping fully (Verde et al. 2002), effectively shut- 
ting off the formation of stars in the first place. 
Gas-rich dwarfs would undergo substantial 
morphological evolution through ram strip- 
ping from the IGM ; their long-term fate 
(destruction or survival) must account for 
such evolution since it will change the dwarf 
binding energy through dissipation (Mayer 
et al. 2001). 

Dynamical friction can in principle provide an 
alternative solution if the clumps spiral in rapidly 
and lose mass owing to tidal heating (Syer & 
White 1998; Tormen et al. 1998). Computer sim- 
ulations and analysis of decaying satellites show 
that a heavy satellite loses up to 90% of its mass 
in a few orbital periods (e.g. Klessen & Kroupa 
1998; Peharrubia et al. 2002; Taffoni et al. 2003; 
see also Hashimoto et al. 2003). van Kampen 
(2002) has argued that the effect of dynamical 
friction may yet be underestimated in computer 
simulations due to limited mass resolution. Dy- 
namical friction, in conjunction with tidal forces, 
will cause the disruption of a satellite after a pe- 
riod of time (e.g. Ibata et al. 1994; Klessen & 
Kroupa 1998; Bullock et al. 2001). Bullock et al. 
(2001) argue that the halo stellar population of 
the Milky Way may be accounted for if sufficient 
satellite dwarf galaxies had already accreted their 
mass at the time the galactic halo formed and are 
then stripped of their less-bound stars by galactic 
tides. Clearly the link between halo morphology, 
halo sub-structure statistics and stellar popula- 
tions offers more avenues for future work. 
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